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We first consider the Einstein-aether theory with a gravitational coupling and a Lagrange mul- 
tiplier field, and then consider the non-minimally coupled quintessence field theory with Lagrange 
multiplier field. We study the influence of the Lagrange multiplier field on these models. We show 
that the energy density evolution of the Einstein-aether field and the quintessence field are signif- 
icantly modified. The energy density of the Einstein-aether is nearly a constant during the entire 
history of the Universe. The energy density of the quintessence field can also be kept nearly constant 
in the matter dominated Universe, or even exhibit a phantom-like behavior for some models. This 
suggests a possible dynamical origin of the cosmological constant or dark energy. Further more, for 
the canonical quintessence in the absence of gravitational coupling, we find that the quintessence 
scalar field can play the role of cold dark matter with the introduction of a Lagrange multiplier field. 
We conclude that the Lagrange multiplier field could play a very interesting and important role in 
the construction of cosmological models. 

PACS numbers: 98.80.Cq, 98.65.Dx 



I. INTRODUCTION 



Ever since the discovery of cosmic acceleration in 1998 
P, , the dark energy has remained a fundamental mys- 
tery, both for its unexpectedly small but non-zero value, 
and for the apparent coincidence of its present density be- 
ing approximately that of the other components. Many 
attempts have been made to address these problems, in- 
cluding, e.g., the cosmological constant, quintessence Q, 
phantom H, quintom [j| , K-essence [f| , holographic dark 
energy M, p , agegraphic dark ener gy M. [l jj], modified 
gravity |13fT Einstein-aether theory |iu7l2l|7and so on. 
It is however fair to say that none of these interesting 
ideas has emerged to be the clear front runner. Many 
of these are only toy models which needs to be further 
developed to be decisively tested by observations [13] ■ 

In this paper, we start our investigation with the 
Einstein-aether field theory [Hill- The Einstein-aether 
theory is an extension of the general relativity theory 
which incorporates a dynamical unit timelike vector field 
A^. The presence of this field breaks the local Lorentz 
symmetry down to a three dimensional rotation sub- 
group. Direct coupling of the aether to the matter would 
violate local Lorentz symmetry yet preserve diffeomor- 
phism invariance. This paper assumes that the aether 
does not couple directly to the matter. 

The Lagrangian density of the Einstein-aether theory 
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is given by [ll| 

% = ^g~[K + \( gfiU A»A» 

where 
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K *I3 = 
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c l9 ^g ia + C2fi?*2 + c 3 5^ . (2) 



K = K a P 1G V a A^pA 



Here Cj and m are constants, g is the trace of the met- 
ric tensor g M „. We emphasize that A is not a constant, 
but a Lagrange multiplier field or auxiliary field. The 
Einstein-aether is similar to the vector tensor gravity the- 
ories studied by Will and Nordvedt [HI]. However, due 
to the presence of the Lagrange multiplier field, there is a 
crucial difference: the vector field is constrained to have 
unit norm. 

Carroll et al. [lH [ItJ showed that the models with 
the generic kinetic term given above are plagued by ei- 
ther ghosts or tachyon, therefore physically unacceptable. 
They found that only the timelike Sigma- model 0, [13] : 



+\( gitu A»A u + m 2 )] , 



(3) 



is well-defined and stable. In this paper, we shall mod- 
ify the fixed-norm constraint from the unit norm to an 
environment-dependent norm. In detail, we modify the 
constraint condition from 



to 



g^A^A" + m 2 = , 



G„ V A»A V + m 2 = , 



(4) 
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where G p „ = R^u — h Rg^v is the Einstein tensor. We will 
shortly find that this modification is not trivial. In fact, 



2 



in our model the energy density of the aether is nearly a 
constant during the entire history of the Universe, while 
the energy density of the unit norm aether is proportional 
to H 2 (H is the Hubble parameter) [l8| . 

Inspired by this interesting result, we also consider 
the Lagrange multiplier field in the non-minimally cou- 
pled quintessence theory. After we completed the discus- 
sions on the non-minimally coupled quintessence case, 
we learned that the Lagrange multiplier field had been 
introduced earlier by Mukhanov and Brandcnbcrger in 
flpj ] . We find that the energy density of the quintessence 
can also be kept at nearly constant value in the mat- 
ter dominated Universe. This suggests a way of finding 
a dynamical origin of the cosmological constant. Fur- 
thermore, by exploring a minimally coupled quintessence 
field with the Lagrange multiplier field, we find that the 
quintessence behaves as pressureless matter, which en- 
ables it to be also a promising candidate of the cold dark 
matter. It is also found that the adiabatic sound speed 
and the rest-frame sound speed of the quintessence is 
exactly zero. In this respect, our conclusion is consis- 
tent with the investigation of Ref. [20(. Different from 
Ref. [20| . we also investigated the Lagrange multiplier 
field in the Einstein-aether cosmology with the modified 
constraint conditions. 

We note that many new results have appeared ever 
since the first version of this work is present in the 
preprint form j2lT[29| . Capozziello et al. [2l[ studied 
the scalar-tensor theory, k-essence and modified gravity 
with Lagrange multiplier constraint. They conclude that 
the well-known mathematical equivalence between scalar 
theory and f(R) gravity is broken due to the presence 
of constraint. Then the important conclusion leads us 
to look for the viable gravity theory among vast origi- 
nally non-realistic ones. Using the scalar and vector La- 
grange multiplier method, Nojiri et al. [22j proposed a 
class of covariant gravity theories which have nice ultra- 
violet behavior and are potentially renormalizable. Cai 
and Saridakis [231 ] investigated the cyclic and singularity- 
free cosmological solutions using the Lagrange multiplier 
method. They showed that the realization of cyclicity 
and the avoidance of singularities is very straightforward 
in this scenario. Feng and Li [24| calculated the primor- 
dial curvature perturbation for the curvaton model in 
the presence of a Lagrange multiplier field. On the other 
hand, Kluson [25| developed the Hamiltonian formalism 
for the Lagrange multiplier modified gravity. For more 
great detail, we prefer the reader to the nice review paper 
Ref. [13|. 

The paper is organized as follows. In section II, we will 
investigate the cosmic evolution and observational con- 
straints on the Einstein- aether with the modified norm- 
fixing condition. In section III, we investigate the cosmic 
evolution of the non-minimally coupled quintessence in 
the presence of a Lagrange multiplier field. In section VI, 
we show the minimally coupled quintessence behaves as 
pressureless matter in the presence of a Lagrange multi- 
plier field, and calculate both the adiabatic sound speed 



and the rest frame sound speed. Section VI discusses 
the results and concludes. We shall use the system of 
units with G = c = h = k = l and the metric signature 
(— , +, +, +) throughout the paper. 



II. COSMOLOGICAL CONSTANT FROM 
EINSTEIN- AETHER 

A. equation of motion 

We start with the following form of action: 



(6) 



The equation of motion obtained by varying the action 
with respect to A enforces the norm-fixing constraint 



G ilv A i A v + m 2 = . 



(7) 



The equation of motion is obtained by varying the action 
with respect to A^, 



V„ {V a A a ) = XG^ a A a . 



(8) 



Multiplying Eq. © with A 1 on both sides, and use the 
constraint Eq. ([7]), we obtain 



A 



rA"V M (V a A a ) 



(9) 



The energy momentum tensor is obtained by the action 
varying with respect to g^ v 30]: 



1 tin ~ z 



2 AW 1 



(V a A a ) 2 +g^A a V a (VpA?) 
2A {i y y) (V a A a )] + A [2V a V (A1 (A v) A a ) 
V a Vp (A a A /3 ) - V 2 (A^Av) 
±A a R a{ll A v) + R^A a A a + RA^Ay 
V^V, (A a A a )+g^V 2 (A a A a ) 



+g»u (G^A'A" 



2 )] 



(10) 



We note that G^ u vanishes in the Minkowski spacetimc. 
It is then follows that m = from Eq. (J7J| and the model 
is reduced to 

2 = >/=5(VpA") a , (11) 
which has been studied in Ref. [30]. 

B. cosmological solution 

Observations reveal that the spatial geometry of the 
Universe is almost flat, so below we consider the spatially 
flat Friedmann-Robertson- Walker (FRW) Universe: 



ds 2 



-dt + a (tf (dr 2 + r 2 dn 2 ) 



(12) 
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with a(t) the scale factor. For such a metric, the vector 
must respect spatial isotropy, at least at the background 
level. Thus the only non- vanishing component of the vec- 
tor should be the timelike component. Using the norm- 
fixing constraint Eq. Q, the components of the vector 
field are simply 



.4" 



= (- 



0, 0, 



(13) 



where H is the Hubble parameter. From Eq. © we then 
find the Lagrange multiplier field A must satisfy 



A = -Li 
3Hdt 




(14) 



where the dot above represents the derivative with re- 
spect to cosmic time t. Using the above two equations, 
we find from Eq. (|10p that the energy density and pres- 
sure of the vector field (and the Lagrange multiplier field) 
are given by 



PA 



in 

T 



and 



PA 




(15) 



#0 ' (i6) 



respectively. It can be easily shown that the energy den- 
sity and pressure satisfy the conservation equation 



PA 



-m( PA + p A ) = o 



(17) 



which is just an expression of the conservation law for 
energy momentum tensor 



rpA;fi _ n 



(18) 



On the other hand, the Einstein equations tell us that 

(19) 



- = --(1 
H 2 2 [ 



w) 



Here w = p/p and p, p represent the total energy density 
and total pressure of the comic fluids ( including aether 
field ). Then the energy density can also be written as 



PA 



m 

T 



~] 2 



(1 



(20) 



For the radiation dominated Universe, we have w ~ 1/3. 
For the matter dominated Universe, we have w ~ 0. 
Hence we have pa — 8.33 m 2 and pa — 6.75 m 2 re- 
spectively. In both cases, the energy density remains 
nearly constant. For the present-day Universe, we have 
w ~ —0.75 such that pa — 3.79 m 2 . In all of these cases, 



the energy density of the aether field behaves as a con- 
stant during the whole history of the Universe. It fol- 
lows immediately that the equation of state of the aether 
is wa — —1. To show the point in great detail, let's 
numerically investigate the Friedmann equation in the 
following. 

Taking into account radiation and matters (including 
ordinary matter and dark matter), we obtain the Fried- 
mann equation 



3H 2 



PrQ 



PmO 
n3 



Til 

~3~ 



(21) 



Here p r o and p m o represent the energy densities of radia- 
tion and matter, respectively, in the present-day universe. 
Define 

x = In a . (22) 
The Friedmann equation can be rewritten as 



3H Z = 8tt 



Let 



p r0 e 



Pmoe 



-3x 



3 I 3 H dx) 



, _ H _ p r0 _ p rnQ 2 _ 87rm 2 

Ho Po Po vH 2 



(23) 



(24) 



then we have 



1 dh\ 



*' = lW-+<Wr" + C (3---I • (25) 

Here Hq, po are the Hubble parameter, total energy den- 
sity for the present-day Universe. £ is a dimensionless 
free parameter. The standard cosmological model, e.g. 
as in Komatsu et al. |3l| . predicts the present matter 
density ratio fl m o = 0.25 and radiation density ratio 
il r o = 8.1 • 10~ 5 . Using this result and taking ( = 0.25, 
we plot the dimensionless Hubble parameter h via red- 
shifts z in Fig. 1 for our model and the standard ACDM 
model. We conclude from the figure that it is consistent 
with the observations. 

In order to show the aether model can lead to the cos- 
mic acceleration, we plot the evolution of deceleration 
parameter q: 



I + 



dlnh 



dx 



(26) 



in Fig. 2. We find the model predicts nearly the same 
transition redshift zt = 0.8 for the Universe from decel- 
eration to acceleration as the standard ACDM model. 

In Fig. 3, we plot the dimensionless energy density 
(from Eq. J25J) 



D = ( 2 



3-- — 

h dx 



(27) 
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FIG. 1: The dimensionless Hubble parameter /i with redshift 
z. The circled line is for the standard ACDM model. The 
solid line is for the aether model. 
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FIG. 2: The deceleration parameter q with redshift z. The 
circled line is for the standard ACDM model. The solid line 
is for the aether model. 



for the aether field with the redshift. It shows the energy 
density of aether field is nearly a constant at redshifts 
greater than 4. 

Carroll et al 18] showed that the energy density of 
aether field is proportional to H 2 in the absence of grav- 
itational coupling. The reason for this difference could 



be understood as follows: for the usual Einstein-aether 
theory, the norm constraint condition Eq. Q gives (by 
order) in flat FRW Universe [3 



A" 



m . 



(28) 



The norm of the vector is fixed as a constant. On the 
other hand, dimensional analysis suggests that we must 
have 



H 



(29) 



Therefore, we conclude that the energy density is 

V tl A ll V u A v ~ H 2 , (30) 

which is consistent with the calculation of Ref. [l8| . For 
our Einstein- aether field, from the constraint equation 
Eq. 



ALL 771 m 



(31) 



where p is the total energy of the Universe. It is apparent 
that the norm of our vector is environment dependant. 
The energy density turns out to be 



H 2 ■ H~ 2 



const. 



(32) 



As another example, we have checked with detailed cal- 
culation that for the Einstein-aether theory with the La- 
grangian density 



Jzf 



-g [G^V^VM Q + A (G m G^A v + m 2 )} 



the energy density of aether in this model is also kept at 
nearly constant level. However, the corresponding Ein- 
stein equations involves derivatives higher than second 
order. 



III. COSMOLOGICAL CONSTANT AND 
PHANTOM FROM QUINTESSENCE 




FIG. 3: The dimensionless energy density D with redshift z. 
At redshifts greater than 4, the energy density of aether field 
is nearly a constant. 



A. equation of motion 

Inspired by the dynamical behavior of the Einstein- 
aether theory with the Lagrange multiplier field, we con- 
sider its role in quintessence field theory. For our pur- 
pose, here we consider the following model: 



-fJ 



, 1 m" 



(34) 



where m is a constant, R is the Ricci scalar. The case 
where A is a constant has been investigated extensively. 
However, we emphasize that A is not a constant, but a 
Lagrange multiplier or auxiliary field. This brings inter- 
esting new possibilities. 
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The equation of motion for A and <f> are given by 

R<p + m 2 = 
A 



(35) 
(36) 



G^VVV - -R = , 



and the energy-momentum tensor takes the form [32-3 

-V^V^D0-i(V<^) 2 G M!/ 
+3m, [ - ^V Q V^ V Q V^ + \{Udpf 
-V«0 V/^JZ *] - \Xm 2 g^ + XG^cj) 
+A 5aw V 2 0- AV^V^ . (37) 



The reason for this difference could be understood 
as follows. The fixed-norm constraint Eq. ([3"5"]) (also 
Eq. ([381)) tells us 



m 
IP 



m 
P 



(44) 



In other words, the strength of the field is determined by 
the energy density of the background matter source. In 
the usual quintessence theory, the strength of the field 
can be arbitrarily large. However, we have ~ H and 
G u 



H 2 



H 2 -H 2 ■ H~ 4 



const. 



(45) 



Another interesting model is 



B. cosmological solution 



-9 



1 / 1 m 2 

-V M 0V^ + A(-i?0+ — 



,(46) 



With the flat FRW background model, we obtain 



nr 



R 6H+12H 2 ' 



A = —(, 



3H 2 



3H12H + ZH' 



(38) 
(39) 



From the energy-momentum tensor, the energy density 
of scalar is 

Pct , = ^H 2 4> 2 + 3XH 2 <j) + 3XH<}> + ^Xm 2 (40) 

Using Eq. ([55)1 and Eq. (|3"9l , and assuming the total equa- 
tion of state w is a constant, we find the energy density 
is 



P0 = 18m 4 • 



(1 + w) (l — w — w 2 ) 
(1 - 3w) 3 



(41) 



We have w ~ for the matter dominated Universe. The 
corresponding energy density is 



= 18m 4 



(42) 



This density does not vary as the Universe expands, so 
the field behaves as a cosmological constant. 

For the radiation dominated case, we have n ~ 1/2. It 
is apparent that Eq. (14"TT) is divergent. This is not surpris- 
ing, since R = in radiation dominated Universe. We 
should look for the energy density from the Lagrangian. 
Ref. [35] showed that the energy density scales as 



Pa 

„2 



(43) 



in the radiation dominated epoch, with po an integra- 
tion constant. Then the field has the equation of state 
w = —1/3. The energy density evolution given in 
Eq. (|4"Tj) is different from the usual non-minimally cou- 
pled quintessence field. For the present-day Universe, we 
have w ~ —0.75 such that ~ 0.16 m 2 . 



in this case the scalar field behaves as "phantom" , with 
equation of state w < — 1. This can be seen by the fol- 
lowing reasoning: the norm-fixing constraint yields 



m 
IP 



So the energy density is 



H 2 ■ H~ 4 



H' 



(47) 



(48) 



The energy density increases with the expansion of the 
Universe, i.e. it behaves as a phantom. By using the 
Lagrange Multiplier field, one may construct the dark 
energy model which exhibit the desired behaviors. 



IV. COLD DARK MATTER FROM 
QUINTESSENCE 

A. the theory 

In Section II and Section III, we have shown that 
the energy density of the Einstein-aether and the 
quintessence can be kept at nearly a constant value in 
the evolution of the Universe by employing a Lagrange 
multiplier field, which points to a dynamical origin of cos- 
mological constant. In this section, we shall show that 
with proper choice of the Lagrange multiplier field, the 
cold dark matter behavior can also be obtained. 

Note that the usual quintessence field can have zero 
pressure for particular potential forms. However, the 
sound speed for the canonical quintessence field is 1, so 
it is not a true cold dark matter. Here we show how with 
the help of Lagrange multiplier field, we can construct 
true cold dark matter with the quintessence field. 

We consider the following Lagrangian density: 



j*f = -Xi/=g 



~V^V> + V (0) 



(49) 
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where A is a Lagrange multiplier field, and V (4>) is the 
scalar potential. The equation of motion for A and cj> are 
given by 



-v^v^ + y(0) = o, 

V M (AW) - XV' = , 



(50) 
(51) 



Here prime denotes the derivative with respect to <\>. The 
energy-momentum tensor is 



Tuv — A 



(52) 



For a flat FRW model, Eq. (|50J and Eq. ((STJl reduce to 
1 ■ 



~<bf -V = 



(53) 



and 



4> + 3H<f) + V' + A _1 A0 = . (54) 
Differentiating Eq. (1531 with respect to time t, we obtain 



V = 



Substituting it into Eq. (fM]). we obtain 



dcj) 3 da 1 dX 
T + + 2~T ~~ 



Solving the equation, we have 

Pa 



A = 



a 3 <fi 2 



(55) 



(56) 



(57) 



with po an integration constant. Using Eq. (|5"2")l and 
Eq. (|57[) . we find that the energy density evolves as 



Pa 



(58) 



which is exactly the behavior of the cold dark matter. We 
note that this expression of energy density is independent 
of the explicit form of the scalar potential. 



B. speed of sound 

A dark matter model will fit the CMB data only if the 
rest-frame sound speed is indeed zero [36J . As this is such 
an important condition, it is worthwhile to spend some 
time to derive the sound speed from the perturbation 
equations directly. For this purpose, we work in the New- 
tonian gauge. In the absence of anisotropic stress, and 
for scalar perturbations, the perturbed flat FRW metric 
can be written in the form 



ds 2 



(1 + 2$) dt 2 + a (t) 2 (1 - 2$) dx i dx i , (59) 



where $ is the Newtonian potential. For the scalar field 
and the Lagrange multiplier field, we define the pertur- 
bation as 



4>{t,x) = 0o (t) + 54>(t,x) , 
X(t,x) = X {t) + SX(t,x) . 



(60) 
(61) 



The perturbed energy-momentum tensor is 

ST° = 5 Plj} = \ (fast - $0 O 2 + V54>) + p^oSX/X , 
ST} = -Spt (i, x) 5\ = -X (foSf - $^o 2 - VSc^j 5\ , 
ikST^ = ik (p0 O + P<po) Sv^ (t, x) = p^V . (62) 

It is well-known that both the adiabatic sound speed and 
the rest frame sound speed (the sound speed for the fluid 
in its rest frame) play important roles in the discussion 
of structure formation theory. Here we work out both 
of these quantities ex plic itly. The adiabatic sound speed 
squared is defined as [33, [38| 



n 2 _ P<P0 
"a — • 

P<j>Q 



(63) 



Since the pressure of this quintessence is vanishing, we 
conclude that 



cl =0. 



(64) 



The rest frame sound speed squared c 2 of the scalar 
is related to the pressure perturbation in the Newtonian 
gauge through [37], HI] 



5 P4> = c? s 5 P4> + -p- (c 2 - c 2 a ) P4>0 y . (65) 



Taking into account of c 2 = 0, we can rewrite Eq. (f6"5)l 
as 



(66) 



In order to calculate Sp^, we resort to the equation of 
motion for the scalar field perturbations. To this end, 
we insert the expansions Eq. (|60H61[) into the Lagrangian 
density, and expand to quadratic order 



(67) 



where «Sfi and ££2 are the Lagrangian terms which are of 
linear and quadratic order in 6(j>, respectively. Variation 
of _Sfo give us the equation of motions for the unperturbed 
(j) and A. Variation of ££2 give us the equation of motion 
for 6<j> and SX, respectively. Variation of does not 
yield new information. 
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Straightforward calculation yields 
_gf 2 = ~a 3 S\(<I>(j)l- 00(50 + V5(j) 
-2a 3 $(5A ( -\^> 2 + V 



-2a^X 



2^<j> l 



2a 



2 (V<ty)' 



(68) 



Variation of ^2 with respect to (5 A and taking into ac- 
count of Eq. (fS"3|) . we obtain 



V 







(69) 



Comparing this with the second equation in Eqs. (|62j) . it 
follows that 



,577 = -6pt(t > 3!)8i=0. 
We then conclude immediately that 
c? =0 . 



(70) 



(71) 



Therefore, for this model the pressure, the adiabatic 
sound speed, and the rest frame sound speed all vanishes, 
which ensures that the scalar field has every desirable 
property of the cold dark matter. 



V. CONCLUSION AND DISCUSSION 

In conclusion, we have introduced gravitational cou- 
pling to the usual Einstein- aether theory and Lagrange 



multiplier field to the scalar held quintessence theory. 
We hnd that by changing the norm-fixing condition from 
unit norm in the usual Einstein aether theory to an 
environment-dependent one, one can make the aether 
density stay constant as the Universe expands. 



As for the case of scalar held quintessence, by intro- 
ducing the Lagrange multiplier field, the strength of the 
quintessence field can be related to the background mat- 
ter source, and the density can also be kept as nearly con- 
stant during the cosmic expansion. This suggests a pos- 
sible dynamical origin of cosmological constant. Alterna- 
tively, we also hnd an example in which the quintessence 
with Lagrange multiplier field behaves as phantom. 



Furthermore, with the help of a Lagrange multiplier 
field, we also proposed a way to generate cold dark mat- 
ter from quintessence. It is found that the pressure, the 
adiabatic sound speed and the rest-frame sound speed of 
the quintessence is exactly zero. These properties enable 
the scalar held to be a potential candidate for cold dark 
matter. 
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